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Abstract 

We present a new model to explain how particles (solar energetic particles; 
SEPs), accelerated at a reconnection site that is not magnetically connected to 
the Earth, could eventually propagate along the wcll-connectcd open flux tube. 
Our model is based on the results of a low-/3 resistive magnctohydrodynamics 
simulation of a three-dimensional line-tied and initially current-free bipolc, that 
is embedded in a non-uniform open potential field. The topology of this configu- 
ration is that of an asymmetric coronal null-point, with a closed fan surface and 
an open outer spine. When driven by slow photospheric shearing motions, field 
lines, initially fully anchored below the fan dome, reconnect at the null point, 
and jump to the open magnetic domain. This is the standard interchange mode 
as sketched and calculated in 2D. The key result in 3D is that, reconnected 
open field lines located in the vicinity of the outer spine, keep reconnecting 
continuously, across an open quasi-separatrix layer, as previously identified for 
non-open-null-point reconnection. The apparent slipping motion of these field 
lines leads to form an extended narrow magnetic flux tube at high altitude. 
Because of the slip-running reconnection, we conjecture that if energetic particles 
would be traveling through, or be accelerated inside, the diffusion region, they 
would be successively injected along continuously reconnecting field lines that 
are connected farther and farther from the spine. At the scale of the full Sun, 
owing to the super-radial expansion of field lines below 3 R©, such energetic 
particles could easily be injected in field lines slipping over significant distances, 
and could eventually reach the distant flux tube that is well-connected to the 
Earth. 
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1. Introduction 



The coronal magnetic field is the main driver of the eruptive solar activity and 
ensures the coupling between the corona and the interplanetary medium. One 
of the manifestation of the solar activity is the acceleration of solar energetic 
particles (SEPs), electrons and ions, at energy ranges from few keV to few GeV 
and their release in the heliosphere. In the Sun's corona, these energetic particles 
can be accelerated either at low altitudes inside magnetic reconnection regions, 
i.e. in flare locations, or at higher altitudes at the shock driven by coronal mass 
ejections (CMEs). 

It is generally believed that cross-field diffusion is weak during the interplan- 



etary transport of energetic particles (Droge et ai, 20101. According to that. 



energetic particles can only reach the Earth if they are eventually injected 
along the very field lines connecting the acceleration site and the Earth. Here- 
after, we refer to the latter as "Earth-connected field lines" . The Solar Particles 
Events are divided in two categories that depends on the time profile of ener- 
getic particle flux measured at Earth : the gradual-SEP and impulsive- SEP. 
The gradual events are usually considered as being accelerated through the 
CME-driven shock, implying a large longitudinal injection of energetic particles 
( [Reames, 1999[). Thus, energeti c particles eventually reach the "Earth-connected 



held lines" ( |Cliver et ai, 2004[ |Nitta et al., 2006D. On the contrary, the impul 



sive events are accelerated through the flare ( Cane, McGuire, and von Rosenvinge, 1986H 



[Reames, 1999[ ) and their interplanetary injection may not be straightforward. In- 
deed, according to several studies ( |Lin, 1970{|Kallenrode, Oliver, and Wibberenz, 1992H 
[Nitta et al, 2006[ ), energetic protons and electrons impacting the Earth may be 
related to flares located far from the footprints of the Earth-connected field 
line. Moreover, it is noteworthy that the flare-accelerated particles are some- 
times detected by different spacecraft in the heliosphere over a wide longitudi- 
nal range around Earth-connected field lines ( |Reinhard and Wibberenz, 1974[ 
[Wibberenz and Oane, 2006] ). 

These fiare-accelerated particles are believed to be accelerated at low altitude 
within the corona and not at the CME-driven shock. Thereby, their injection into 
the interplanetary medium should be achieved through a reconnection process 
between open and closed coronal magnetic field. This reconnection can either be 
the same as, or different from, the flare reconnection. 

This paper addresses the following questions : How can magnetic reconnection 
allow energetic particles that are accelerated in the low solar corona, far from 
the "Earth-connected field lines" , to get access to the Earth ? Do exist dynamics 
of the magnetic reconnection that would lead to inject fiare-accelerated particles 
over a wide angular distance in the heliosphere ? Hereafter, we link the three 
issues of the open versus closed corona, of particle injection problem and of 
three-dimensional reconnection, in order to motivate the model that we propose 
to answer these two questions. 

1.1. The Coupling Between the Closed Corona and the Heliosphere 

The large-scale coronal magnetic field presents two different types of connec- 
tivity: closed and open. In the closed magnetic domain, field lines have both 
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footpoints rooted in the photosphere, while for the open domain, a field line 
is anchored at the photosphere by one footpoint only, its other end extending 
into the heliosphere. The coupling between the corona and the interplanetary 
medium is carried out at the transition between open and closed magnetic field. 
This coupling can either be ideal {e.g. when propagating waves carry energy 
from one domain to another) or non-ideal {e.g. when magnetic reconnection 
leads to mass and magnetic flux transfer between both domains). Examples 
of coronal magnetic configurations displaying such an open/closed transition 
are : at large-scales, the helmet streamers and plumes/pseudo-strcamers {e.g. 



Wang, Sheeley, and Rich 120071 Wang et al. I2007P : at smaller scales, coronal 
lets (e.o. lCirtain et aLll2007p . Streamers are formed above bipolar regions, where 
the closed coronal loop system is opened at high altitudes to the interplanetary 
medium through the solar wind pressure {e.g. [Grappin, Leorat, and Habbal[ 
I2002p . Pseudo-streamers, and jets respectively, are created when a larger active 
region {see e.o. lDel Zanna et fflZ.[l201I1 ) or a small magnetic bipole (Shibata et al., 1994H 



Moreno-Insertis, Galsgaard, and Ugarte-Urra, 2008 Torok et al, 2009 1 emerges] 



in an already open coronal magnetic flux tube. 

It has been shown that, apart from special cases (Titov et al., 20101, a single 
bipole embedded in an open fleld most often relates to a magnetic topology 
that comprises a single magnetic null point in the corona ( as in lAntiochosl 
[TM51 [Pariat, Antiochos, and DeVore| [MM ITorok et al.\ [20®- Such a null- 
point is associated to a fan separatrix surface and a singular spine fleld line 



( Lau and Finn, 1990 1. For a bipole embedded in an open field, the fan forms a 
dome which encloses an inner connectivity domain. The fan separatrix surface 
separates closed magnetic field lines from the surrounding open ones, located in 
the outer connectivity domain. The spine field line crosses the fan surface at the 
null point. Therefore, its outer part is also open. 

Magnetic reconnection at the null point allows the exchange of connectivity 
between closed and open magnetic field lines. This is the standard "interchange 



reconnection mode" ( Pariat, Antiochos, and DeVore, 2009|[Edmondson et al., 2009[ ).[ 



In general, this closed to open reconnection allows a direct coupling between the 
open and the closed corona. In particular, it is a natural mechanism to inject 
impulsive-SEP related energetic particles from the closed corona to the interplan- 
etary medium. The acceleration may occur either in a fiare located inside the 
inner domain, before the null point reconnection itself, or at the null point during 



the interchange reconnection (Litvinenko, 2006 Rosdahl and Galsgaard, 2010 



Browning et al., 2010). The detection at the Earth of such flare-accelerated par- 
ticles implies that the open coronal flux tube surrounding the outer spine and 
involved in the interchange reconnection, must be a priori directly connected to 
the Earth-connected field lines. But, it is not always the case. 

1.2. The Problem of the Particle Injection 

Generally, the geometry of the interplanetary magnetic field (IMF) is assumed to 
follow the Parker spiral. Depending on the solar wind velocity, this implies that 
the Earth-connected field lines are rooted between 40° and 70° in longitude. 
However, past studies of impulsive-SEP events show that energetic particles 
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can very often be detected at the Earth even if the associated flare is located 
away (up to 50°) from the Earth-connected interplanetary spiral field lines 
dKallenrode, Oliver, and Wibberenz, 1992[[Nitta et al, 2006[ ). This may first im-| 
ply that the Earth-connected field line is not a simple Parker spiral starting in 
the low corona. 

Indeed, using PFSS extrapolations iKlein et aZT](|2008p showed that open mag- 
netic flux tubes in flaring active regions associated to SEP events, can spread 
in longitude from 4° at the photosphere to more than 70° at the solar source- 
surface, at the base of the Parker spiral. Comparing coronal magnetic field, 
metric radio observations and energetic particles in-situ measurements, they 
showed that the large longitudinal extensions of open magnetic fiux tubes could 
ensure a direct path to Earth for energetic particles, being accelerated far from 
the Earth-connected field lines. Even though this effect is surely important, it 
may not be the only way to inject energetic particles onto the Earth-connected 
field lines. 

In addition to this first issue, energetic particles of some impulsive SEP events 
can also be detected in a wide flux tube, over more than 90° of longitude 
from the the parent active region ( [Wibberenz and Cane, 20"06l ). These large 
particle beams can be explained by lateral transport along coronal magnetic 
loops (Newkirk and Wentzel, 19781 or lateral diffusion in open coronal magnetic 
field dReid, 1964 Wibberenz et ai, 1989[ ). Instead invoking the transport effect 
of particles, could this large angular spread result from a specific dynamics of 
injection channel related to the interchange reconnection mode? 

For conventional models of interchange reconnection, we expect that parti- 
cles can only be injected into a narrow bundle of field lines whose connection 
to the Earth is rather unlikely. Indeed, the field line that reconnects at the 
null point and opens in the outer connectivity domain is located where the 
outer spine was located at the time preceding the reconnection. Thus, as the 
reconnection evolves, the reconnected field lines form an open flux tube concen- 
trated around the initial location of the outer spine. Flare-accelerated particles 
should therefore be only injected into this localized flux tube. Even though 
the open magnetic field can display a rapid (super-radial) divergence in the 
corona (as easily seen in |Wang, Sheeley, and Rich[ 120071 IKlein et'oH I2008p . 
particles are transferred into the open corona along this localized flux tube, 
implying a narrow localized particle beam. Nevertheless, the spine itself can 
shift its position by several degrees during a long-enough reconnection process 
(see eg. [Moreno-Insertis, Galsgaard, and Ugarte-Urra[ [2008; Ednionds on et all 
I2009p . and could lead to broaden the open flux tube along which particles 
would be injected. But, we ignore if this motion is directed towards the Earth- 
connected field line and if the displacement covers a sufficient angular distance. 
This connectivity problem is illustrated by the theoretical geometry sketched 
in Figure [T] This paper addresses the following question : how the dynamics of 
magnetic reconnection in an open-spine null-point topology can lead to inject 
energetic particles in the interplanetary medium over a large longitudinal range? 
Through the invoked mechanism, the angular spreading can lead particles to 
eventually reach the Earth-connected field line, even if the outer spine will never 
be connected to it. 
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1.3. The Role of the Magnetic Reconnectioii Dynamics 

One first clue to answer this question can be obtained from the analysis of a for- 
mer 3D interchange reconnection model driven by flux emergence ( )Rosdahl and Galsgaard, 2010[ ) 
Particle acceleration and transport was there modeled by following test-particles 
in a fixed magnetic field taken from one snapshot of an MHD simulation. They 
showed that particles are accelerated at the reconnection current sheet, and are 
injected along a long narrow layer of open field lines. The spatial extension of 
this propagation channel could be a way to get energetic particles in a large 
range of longitudes. One can still ask why such an extended acceleration region 
occurs in this simulation, and whether or not it is sufficient to explain the 
large spreading of energetic particle beams. According to lAulanier et aL\ (I2007p 
and lMasson et aL\^'2009\\ the injection channels arc affected by the dynamics of 
magnetic reconnection. Therefore, the time evolution of the magnetic field could 
play an important role in the extension of particle beams. 

So as to address all these questions, we performed a new 3D MHD simulation 
of a coronal null-point topology with an open outer spine, driven by slow pho- 
tosphcric shearing motions, and we investigate the dynamics of the interchange 
reconnection. The model is not dimcnsionalizcd, but the topology being consid- 
ered can make it relevant to a pseudo-streamer around an active region, or to a 
small jet inside a coronal hole. After describing the numerical model in Section[51 
we present the results on the magnetic reconnection dynamics in Sections [3] 
and m Assuming that the dynamics of the magnetic field can represent the 
macroscopic evolution of injection and propagation channels, our results provide 
a new conceptual model to understand how energetic particles could reach the 
Earth-connected interplanetary field line when the reconnection site is initially 
located far away from it (Figure [T|). We finally discuss the relevance of our model 
on the interplanetary injection of flare-accelerated particles (Section [Sj . 

2. Model Description 

The initial magnetic conflguration is potential and is built with sub-photospheric 
monopolcs placed such as to obtain a bipolc embedded in an open diverging 
magnetic field, leadinan asymmetric null-point topology with an open outer 
spine. The configuration is therefore well-suited to model large-scale pseudo- 
streamers above active regions, or small-scale jets inside coronal holes. The 
magnetic field is normalized so that the maximum initial photospheric magnetic 
field B[z = 0,t = 0] = 1. The top left panel on Figure [5] displays the initial 
magnetic topology where grey field lines materialize the null-point separatrices. 

We use the explicit 3D MHD code described in |Aulanier, Demoulin, and Grappinj 
(|2005p . to which we added gravity finite-/? effects. The numerical simulation is 
performed in non-dimensionalized units and in a Cartesian domain in which x 
and y are the horizontal axes and z the vertical axis, covering the volume defined 
by: X G [—37,23], y G [—30,30] and z S [0,100]. The mesh is non-uniform with 
201 X 181 X 201 grid points, not counting two layers of ghost cells around each 
face of the domain. The mesh size results from a combination of geometric series 
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in all directions and is adjusted to be the thinnest around the null point position 
(z = 4.81) and at the inversion line (z = 0). The mesh is displayed on Figure [2] 
: on the top right panel, the black grid represents the cells of the mesh at the 
photospheric boundary, whereas the lower part {z = [0, 12]) of the vertical mesh, 
where the fan-spine separatrices are located, is shown on the bottom panel with 
the white grid. (The details of the mesh can be seen by zooming in the electronic 
version of this figure.) 

Our calculation solves the following set of equations using a fourth-order 
finite-difference and third-order predictor-corrector scheme: 

|^ = -V.(pu) + «;A(p-p(t = 0)), (1) 
= -(u • V)u + p-^jxB- p-^\/P + g + iy' Du, (2) 



du 

dt 



dT 
'dt 



-(u • V)T -(7-l)rV-u-|-K A(r - T{t = 0)), (3) 
— =Vx(uxB) + r;AB, (4) 



J = V X B , F = 2pT (5) 

p being the mass density, T the plasma temperature, u the plasma velocity, 
B the magnetic field. 7 is the ratio of specific heats chosen here at 1. These 
equations are written in the code in their non-conservative fully developed form. 
The lower boundary of the numerical box follow the line-tied conditions, whereas 
the upper boundary and the edges of the box follow the open condition. 

The operators v'Du and TyAB correspond respectively to a viscous filter that 
uses a pseudo-Laplacian adapted to the local mesh ( [Aulanier, Demoulin, and Grappin, 2005D , 
and a standard coUisional uniform resistive term. The values of v' and 77 lead to 
diffusive speeds such as: = 7% c™^'^ and = 5% d^^^ on the scale of the 
smallest mesh, where the maximum Alfven speed is c'^^^ — 7. k A(p — p{t ~ 0)) 
and kA{T — T(t = 0)) are non-physical explicit diffusive terms which smooth 
the gradients of density and temperature at a speed of = 7.5% c™^'^, and help 
to stabilize the numerical computation. 

Owing to the high-order spatial and temporal scheme of our 3D MHD code, 
firstly the numerical diffusion is very small, and secondly current sheets al- 
ways remain resolved on more than one grid point. Indeed, the use of the 
explicit resistive term is mandatory to stabilize the code (see the Appendix 
of [Aulanier, Pariat, and Demoulin] (j2005p ). Also, the use of five grid points to 
calculate spatial derivatives imposes that all the gradients developing throughout 
the simulation are resolved on no less than three grid points (see the bottom 
panel of Figure [5]). The explicit resistivity in our code is therefore responsible 
for all non-ideal processes described in this paper. 
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In order to drive magnetic reconnection, we prescribed a slow photospheric 
flow which moves the positive polarity toward the negative y in the area of 
X £ [—7,3] and y g [10, —30], where a part of the fan is rooted. The top right 
panel on Figure [5] displays the location at the photospheric boundary and the 
shape of the applied shearing motion. A temporal ramp is applied, allowing an 
initial relaxation phase between t = and t = 15. The photospheric driving is 
then constant until i = 40, after which it is progressively stopped with a similar 
temporal ramp. The maximal amplitude of the prescribed photospheric velocity 
is: u^^^^ = 0.35 w 5% c^'^'^. 

Considering the large size of the numerical box and the existence of a null 
point, a finite-/? stratified atmosphere is required so as to respect the following 
coronal regime : /3 < 1 (except in a radius of 1 around the null point, where 
/3 > 1) and a sub-Alfvcnic and sub-sonic regime in the whole domain for the 
coronal perturbations, directly induced by the shearing motions. 

We define this atmosphere as follows. We first fix the photospheric density 
p[z = 0, t = 0] in order to obtain a maximal photospheric Alfvcn speed c'^'^^ = 7. 
This characteristic speed defines the time unit required for an Alfven wave to 
travel the distance between the centers of the two polarities of the small bipole, 
D = 7. 

We then define the gravity and temperature profiles such that the Alfven 
speed remains roughly constant above the altitude of the null point. There, 
the density scale height, and therefore the pressure scale height = 2T/g should 
be half of the magnetic scale height = B/{dB/dz). Assuming an hydrostatic 
atmosphere, and choosing 2T/g as shown in Figure [3j we integrate numerically 
the logarithm of the pressure and deduce the vertical density profile at t = 0. 
Even though this atmosphere does not exactly reproduce solar numbers, it still 
puts the system in a solar-like regime. 



3. Null-Point Squashing and Magnetic Reconnection 
3.1. Current Sheet Formation 

The photospheric shearing motion of one part of the fan's footpoints leads to a 
bulging of field lines rooted below the fan surface. This expansion compresses 
the fan separatrix surface and the inner spine, leading to the development of an 
over-density (Figure S]). 

This also leads to the shearing of the spine (bottom row, FigureU]), i.e. a loss of 
alignment between the inner and the outer spine ( [Pariat, Antiochos, and DeVore, 2009[ ) .| 
This deformation of the null-point compresses and extends the region of weak 
B, inducing the formation of a narrow and intense current sheet in the vicinity 



of the null point and along the separatrices ( Low, 1987] Aly, 1990 |Lau, 1993} 



Rickard and Titov, 1996 Galsgaard, Priest, and Titov, 20031, as seen on Fig- 



ure ID 

Thus, null-point reconnection is expected to occur, leading to a realignment 
of the inner and the outer spine ( [Antiochos, Karpen, and DeVore, 2002p . By 
plotting field lines from fixed footpoints and by integrating these field lines up 
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to their conjugate footpoints at different times, we can follow the evolution of 
the connectivities during the simulation. 



3.2. Interchange Reconnection 

In order to follow the field lines that are initially rooted below the fan surface, 
in the positive polarity advected by the applied flow, we computed analytically 
the position of their advected footpoints at each times. Therefore, we follow 
the connectivity of these field lines from their footpoints fixed in the applied 
photospheric fiow. We initially plot three groups of field lines: pink, light blue 
and red, closed below the fan surface, and four groups of field lines (pink, light 
blue, yellow and red) open in the corona (bottom left panel of Figure . 

In each group, field lines are plotted from fixed footpoints distributed along 
a small segment. All these segments have the same length and each of them is 
located at the same distance from the fan's footpoints. The field lines of the 
same color group are plotted from footpoints uniformly distributed along one 
segment. As reconnection proceeds, magnetic fiux is exchanged between con- 
nectivity domains, and the fan surface is displaced and sequentially crosses the 
fixed footpoint of the plotted field lines. Since the plotted field lines are uniformly 
distributed along each segment, the number of field lines that reconnected per 
time unit is proportional to the reconnection rate. 

As the system evolves, the field lines, initially closed in the inner connectivity 
domain, reconnect at the null point and jump to the outer connectivity domain, 
opening- up in the corona (bottom row of Figure On the opposite, the initial 
open field lines reconnect at the null point and closing-down below the fan 
surface. This exchange of connectivity between the inner and the outer connec- 



tivity domain is, in fact, the interchange reconnection (Edmondson et ai, 20091. 
During the reconnection through separatrices, the transfer of magnetic flux be- 
tween the inner and the outer connectivity domain leads to displacements of the 
separatrices. Therefore, the positions of the fan surface and the spine line, both 
change. 

Even though the deformation of separatrices appears clearly during the evo- 
lution (see the dark blue field lines lines on Figure S]) the relative length over 
which the spine field line shifts by reconnection, as compared to the horizontal 
scale over which the neighboring open field lines expand with height, is not very 
large. In itself, this result suggests that the large expansion with height of field 



lines above flaring active region, as identified by (Klein et ai, 2008), may not 
always allow energetic particles previously accelerated low down in the corona 
to reach the Earth-connected field lines. 



4. Slip-Running Reconnection of Open Field Lines 
4.1. Slipping Field Lines 

Before field lines reconnect at the null point, they display an apparent continuous 
slipping motion toward the inner spine (see Figure 0] as well as the animation on 
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the electronic version). Figure [S] also displays the evolution of the connectivities 
at high altitude, in the whole numerical domain. The four panels show the same 
groups of field lines, as plotted in Figure at four different times. The animation 
shows that, right after having reconnected at the null point and jumped to the 
outer connectivity domain, field lines are located all along the outer spine in the 
whole numerical domain. Nevertheless, these reconnected field lines do not stay 
in the vicinity of the outer spine: they very quickly move farther away from it. 
Indeed, after the null point reconnection, field lines display an apparent slipping 
motion along a specific direction towards the negative y. 

Even though the topology is that of a null point, this apparent slipping 
motion is reminiscent of a slipping reconnection regime occurring through quasi- 
separatrix layers ( Aulanier et al., 2006p . A quasi-separatrix layer (QSL) is a 



purely three dimensional geometrical object, formed by narrow volumes of con- 
tinuous connectivity but across which the connectivity displays strong gradi- 
ents ( [Priest and Demoulin, 1995| . Contrary to the null-point reconnection, the 
absence of magnetic discontinuity implies that, in a QSL, field lines simply 
exchange their connectivity with their neighbors and that their footpoints is 
moved in time across the QSL. Thus, continuous reconnection between neigh- 
boring field lines leads to a continuous exchange of connectivity and apparent 
slippage of field lines plotted from fixed footpoints can be noticed. This slipping 
reconnection regime has been confirmed by several MHD numerical simulations 
dPontin et al., 2005} [Aulanier et al., 2006P and supports numerous association 
of observed flare ribbons with QSL calculated in extrapolated magnetic fields 
( [Demoulin et al., 19971 [Mandrini et al., 1997 ). 



Since it has recently been shown that this regime can also occur around 
null points with completely closed magnetic separatrices ( [Masson et al., 2009[ 
[Torok et al., 2009| , we investigate whether or not the same happens here, in the 
presence and around open separatrices. The slipping motion of reconnected field 
lines, before and after the null point reconnection, could result either from slip- 
running/slipping reconnection through QSLs, or from mass motions induced 
by the propagation of shear Alfvcn waves injected cither by the photospheric 
driving, by the reconnection jets, or even from numerical errors. 

Note that at late times, the animation shows that the green field lines, rooted 
in the photospheric driving region far from the null point, also move at the top 
of the numerical domain (Figure. [S|). Nevertheless, these motions are straight- 
forward to interpret. They clearly result from the upward propagation of undu- 
lations, that simply correspond to the propagation of a shearing Alfven wave 
from the line-tied boundary, which is generated by the photospheric driving. 

4.2. Relation Between Plasma and Field Line Velocities 



We compared the slipping velocities of field lines and the plasma velocities in 
the vicinity of these reconnected field lines, at the top of the box. We selected 
the fifth reconnected field line in each of the 3 studied groups of field lines, and 
we determined for these three selected field lines their slipping velocity at the 
top of the box for each time interval of the simulation. Then, we evaluated the 
plasma velocity at the top of these three selected field lines. Since the mesh of 
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the plotting is much tighter than the mesh of the 3D MHD code, the position 
of the top of the field lines is never located exactly on one of the mesh points 
where the MHD variables are integrated. Thus, we computed the mean value of 
the plasma velocity of the four mesh points surrounding the field lines position 
at the top of the box. The resulting velocities are plotted in Figure |6l 

After the field lines have reconnected at the null point, their slipping veloc- 
ity is greater than the Alfven speed, ca — 1.5, and than the fast-mode speed 
cf — 1.7 at the top of the box. This leads to a super-Alfvcnic and super-sonic 
slipping motion. The slipping velocities of field lines later decrease to w 0.2 and 
become sub-Alfvcnic and sub-sonic. Meanwhile, the plasma velocity at the top 
of the box is found to increase during the simulation, in agreement with the 
upward propagation of the shear Alfven wave from the photospheric driving at 
the bottom of the box. Even though the slipping velocities decrease, they remain 
much faster than the plasma velocity until t = 50 (see Figure^]). After t = 50, the 
temporal evolution of the slipping velocities becomes the same as the evolution 
of the plasma velocities and their values becomes comparable. Therefore, the 
displacement of field lines at the end of the simulation is indeed due to mass 
motion, i.e. ideally evolving field lines, such as the green field lines. 

At this stage, the analysis suggests that the slipping motions of the recon- 
nected field lines, at least until t ~ 50, are due to magnetic reconnection through 
a QSL, and at least not to mass motion. The evolution of the slipping veloc- 
ity of the field lines, successively super-Alfvcnic then sub-Alfvenic, corresponds 
to a transition from slip-running to slipping reconnection following null-point 
reconnection, as first put toward by IMasson et al.\ (j2009p , but for fully closed 
magnetic fields. 

4.3. Can Numerical Errors Induce the Slipping Motions ? 

If the slipping motion of the reconnected field lines was due to systematic nu- 
merical errors throughout the whole integration domain, whether coming from 
the MHD computation or from the field line integration code, the green field 
lines, which evolve through the upward propagating shear Alfven wave, should 
display a slipping motion too. The animation attached to Figure [5] clearly shows, 
however, that the reconnected field lines slip long before the top of the undulating 
green lines start to slip, the latter only being due to the Alfven wave passing 
through the top boundary. This difference indicates that there is no significant 
systematic numerical error. 

Still, one expects larger errors which could locally affect the reconnected field 
lines, since those pass through a current sheet, in which sharper magnetic gradi- 
ents develop. Such errors can be estimated with the divergence of the magnetic 
field, since our MHD code does not enforce the solenoidal condition V • B = 0. It 
should be noted that all the MHD equations are fully developed in the code, with 
all the V-B terms being removed, so that these terms neither produce momentum 
nor induction. In order to quantify the effect of non-zero V • B on the slipping 
motions, we nevertheless study the dynamics of the field lines with modified 
magnetic field B*. At each time and position in the domain, every component of 
the modified magnetic field was computed by the following Taylor-like expansion 
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where di is the local size of the mesh interval in the direction i : B* = Bi+di V-B, 
where i = x,y,z. Using this modified magnetic field, we re-conduct the same 
velocity analysis as described above, for the very same reconnected field lines. 
We find that at the top of the box, the temporal evolution of the slipping velocity 
and the displacement of the B* field lines display the same temporal evolution 
than the B field lines. The dynamics of magnetic rcconncction for B"^ is the same 
than B : the succession of null-point, slip-running and slipping rcconncction 
regimes. However, one notices that reconnection process for B* starts twice 
Alfven time later than for B. In addition, the slipping velocity of field lines, 
at a given location at the upper intersection of field line with the boundary, 
differs between B* and B. Indeed, we evaluate that the slipping velocities with 
B* depart from those with B by 35% on average. This number is presumably 
due to the strong sensitivity of the integration of field lines within a strongly 
diverging flux tube such as a quasi-sepatrix layers (which is described below as 
the structure along which field lines are slipping), even when starting from a 
region where B has locally been slightly modified to B* 

Even though our MHD code does not ensure the solenoidal condition, and in 
spite of the fact that the present calculation was ran in single precision, we thus 
find that the measurable numerical errors do not affect the large-scale dynamics 
of magnetic reconnection and the observed slipping of the reconnected field lines. 



4.4. An Elongated QSL around the Open Spine 



In IMasson et al.\ (j2009p . we had shown that the closed separatrices of a null 
point were surrounded by a QSL halo. The presence of a QSL, through which 
continuous magnetic reconnection is known to occur, can be checked by com- 
puting the norm N ( [Priest and Demoulin, 1995] ) or the squashing degree Q 
( [Titov, Hornig, and Demoulin, 2002 1. These quantities somehow measure the 
gradients of field line connectivity, and thus they map the position of the QSL. 

Usually N and Q are calculated in closed field configurations only. Here we 
compute Q in our open fields in a similar way as how IDemoulin and PriestI 
(|1997p calculated N : the upper and bottom boundaries are considered as the 
two reference planes in which field line footpoints are calculated. Even if the lower 
line-tied boundary appears as more natural reference plane than the upper open 
frontier at z ~ 100, Q is uniquely defined for each lines given this choice of 
computation limit. However, the choice of a different reference plane (such as 
a different altitude for the top boundary) would result in a different value of 
Q for a given field line. Nevertheless, assuming a monotonous divergence of the 
magnetic field, the distribution of Q would remain monotonously homogenous 
with a different choice of z. Hence, the occurrence and locations of QSLs in open 
field lines, manifested by open narrow layers of high-Q, would be independent 
of the choice for the reference upper plane. 

Figure [7] shows the location and shape of the high-altitude QSL which sur- 
rounds the open spine, as seen from the large Q isocontours plotted at z = 99.5. 
Initially, the QSL is broad and is almost symmetrically distributed around the 
spine. At later times, however, the prescribed photospheric motions shear the 
system, and increase the asymmetry of the null point. Following IParnell et al.\ 
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P996p . the parameter p which defines tlie asymmetry of the fan, increases in 
our simulation from 1.2 (at t = 0) to 2 (at t = 65). This is associated with 
the development of a non-symmetric, elongated and thin QSL around the open 
spine, in a similar way as what we found earlier in an initially asymmetric, but 
closed, null-point topology ( Masson et al., 2009 ). 

On Figure ini we report the distribution of the high value of Q, and thus the 
location of the QSL, at the top of the box. It is now clear that the intersections 
of the slipping field lines with the upper boundary of the box are located within 
the QSL. This finally leads us to interprete the super- (followed by sub-) Alfvcnic 
field line motions as a substantial amount of slip-running (followed by slipping) 
reconnection, of the field lines which have previously jumped from the inner to 
the outer connectivity domain by the interchange null-point reconnection. 

The distribution of Q (Figure [T]) can also explain why the slipping velocity of 
some reconnected field lines displays a local minimum around t = 40 after the 
slip-running phase (see Figure [6]). A band of low Q has formed within a higher 
Q region. lAulanier et al. \ (|2006p noted that field lines slip slower (res. faster) in 
regions of lower (resp. larger) Q . We also find that the brief period of lower field 
line velocities around t = 40 is indeed associated with the crossing of this local 
low-Q region by the slipping field lines. 



5. Discussion 

5.1. MHD Results 

So as to explain several in-situ measurements in the heliosphere and at Earth, 
we investigated an MHD model for the dynamics of 3D magnetic reconnection to 
address the interplanetary injection of energetic particles during impulsive-SEP 
events. The model is based on the interchange reconnection mode in an open- 
spine null-point topology ( [Edmondson et al., 2009[|Pariat, Antiochos, and DeVore, 2009[ ).[ 
This situation is relevant for the two solar magnetic configurations that couple 
the closed magnetic corona with the open corona and the heliosphere : the 
large-scale pseudo-streamers above active regions, and the small-scale jets inside 
coronal holes. In a pseudo-streamer configuration, flare-related particles may 
have been previously accelerated inside the inner connectivity domain, whereas 
in a jet configuration, particles are supposed to be directly accelerated at the null- 
point. In both cases all energetic particles pass through the current sheet formed 
at the interface between open and closed field regions, before being injected in 
the interplanetary medium through interchange reconnection. 

The main result of the MHD simulations is that, before and after the standard 
interchange reconnection at a coronal null point, the field lines rooted in the 
photosphere continuously reconnect through slip- running and slipping reconnec- 
tion regimes. This is due to the presence of a halo of quasi-separatrix layers 
(QSLs) that surrounds the fan separatrix surface and spine singular field line 
originating at the asymmetric null point. In such an "hybrid topology" magnetic 
reconnection does not behave as usual, i.e. field lines do not exchange their 
connections instantaneously by pairs. This continuous reconnection associates 
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with a large extension of the current sheet along the fan, exceeding the length 
of the current sheet which forms at the shearing null point. This large current 
sheet is a direct result of the development of a QSL around the spine line of 
a null point becoming increasingly asymmetric. This QSL sustains slip-running 
and slipping reconnection, hence leads to a large displacement of slipping field 
lines at high altitudes in the open field domain. This result generalizes to open 
spine configurations the same MHD process originally discussed for a fully closed 
magnetic fields ( [Masson et ai, 2009[ ). Since slipping reconnection across QSLs 
around an asymmetric (also called improper) null point seems to occur for line- 
tied cases as well as for open boundary conditions, we argue that it must be 
a generic property of 3D reconnection at null points. As such, the "torsional 
fan reconnection mode" found in asymmetric nulls by lAl-Hachami and Pontini 
(|2010p may actually be just another formulation of the same physical effect. 

Even though our calculations do not treat the dynamics of the particles them- 
selves, our analyses of the evolution of reconnected field lines can be used as a first 
proxy to determine the time-evolution of injection and propagation channels for 
accelerated particles. Thereby, we assume that the particle acceleration process, 
wherever it occurs, can lead to a significant number of particles at sufficiently 
high energies. 

5.2. Possible Consequences for Particles 

In the classical two-dimensional interchange reconnection mode, particles should 
only be injected into open field lines through the short current sheet that forms 
around the misaligned separatrices of the null point. Thus, particles should only 
be injected along the few field lines which have reconnected at the null point. The 
QSL-related field line dynamics in our simulation, however, let us conjecture that 
during 3D interchange reconnection, particles can also be successively injected 
along (initial fast) slip-running (and later slow slipping) field lines. 

In principle, the association of slip-running and interchange reconnection 
should allow that particles are accelerated all along the intense current sheet 
localized in the QSL and not only around the null point. Let us illustrate a 
discretized version of this continuous mechanism : when an initially closed field 
line (i) reconnects and opens at the null point, a first set of particles are 
injected along this reconnected open field line. Then, as the same open field line 
(i) slip-runs away from the open spine, another closed field line (« + l) reconnects 
and opens at the null, and a second set of particles (i -I- l)i are injected along 
this second reconnected field line. In the meantime, the motion of the field line 
(i) has made it slip along the QSL-related low-altitude current sheet and a new 
set of particles(i)2 is injected along the open field line (i). As this whole process 
repeats itself n times, newer and newer sets of particles are accelerated and 
injected along the field line (« -I- n) being just reconnected at the null point, as 
well as along all the slip-running field lines {i + n — l]i + n — 2; + that are 
moving farther and farther away from the open spine. Therefore, this mechanism 
implies that the set of energetic particles would be injected continuously along 
field lines which are undergoing slip-running/slipping reconnection away from 
the open spine. 
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We conclude that during the interchange foUowed by shp-running recon- 
nection, particles could be injected within an extended and narrow flux tube, 
progressively formed by the field lines slipping along a large-scale open QSL. In 
a fixed magnetic field configuration, this may explain the elongated shape of par- 
ticle impacts at the top of the numerical domain, as obtained in the simulations 
of |Rosdahl and Galsgaard] (|2010p . In a time- varying system, our analysis implies 
that particle beams originating from the interchange reconncction should sweep 
the high altitude domain, along a narrow QSL-related flux tube that extends 
through the slipping reconncction process. 

5.3. Adaptation to the Heliospheric Context 

When considering a solar-like configuration such as a large-scale pseudo-streamer 
or a small-scale jet , the resulting sweeping energetic particle beam could even- 
tually reach the Earth-connected Parker spiral field line, even if the acceleration 
site is not initially connected to this very line (as shown in Figure [1]). Our model 
provides a new type of MHD-based particle injection process along a sweeping 
beam, that may cover a large angular range. 

Around 80 % of active regions related to impulsive SEP events are close to or 
included in an open coronal magnetic fiux tube ( Wang, Pick, and Mason, 2006[ 



[Pick et al., 2006| . They are auspicious magnetic configurations for finding pseudo-| 
streamer and jet configuration. In addition, in such a region nearby open coronal 
magnetic field, when we measured a SEP events at the Earth, a small jet is 
usually observed in UV and occurs at a comparable time than type III bursts 
dWang, Pick, and Mason, 2006| |Nitta et al, 2006D .This suggests that open and 



close field lines are reconnecting through the interchange mode in a null-point 
topology. Therefore, the null-point configuration embedded in an open magnetic 
field, as shown in Figure [TJ can be considered as a typical configuration for flare- 
related impulsive SEP events, regardless any CME. This leads us to suggest 
that the model proposed by this study should be a common injection process. 
However, one can ask if the angular distance travelled by slipping field lines and 
thus the sweeping particle beam, in a solar context, is significant ? 

Even though the non-dimensionalized units do not allow to directly compare 
the results of the simulation with the observations, one can estimate the angular 
distance, 9, covered by slipping field lines for typical solar values. Associating the 
geometry of the system to a right triangle, where the spine and the top of the box 
form the right angle and the slipping field line is the hypothenuse, we evaluate 
the angle 9. For an Alfven speed of ca ~ 500 km s~^, and an arbitrary distance 
unit of = 10 Mm for a coronal jet (resp. 100 Mm for a pseudo-streamer), 
the fiux tube formed by the slipping field lines expands over an angular range 
6* ~ 14 ° at the altitude of 142Mm for a jet (resp. 2 R©). Comparing t the spine 
field line's motion covering the angular distance ^spino — 1-4 °, it appears that 
the slipping reconncction lead to a much more substantial angular spreading of 
the open coronal flux tube, guiding energetic particles in the heliopshere. 

In addition, it is noteworthy that 9 is constrained by the simple geometry 
of the the model used in the 3D-MIID simulation. Indeed, the diverging of the 
magnetic field in 3D is relatively small, but a largest expansion will lead to 
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increase values of 6. Actually, the spatial extension of the open coronal magnetic 
field being super-r adial ( |Wang, Pick, and Mason, 2006| [Klein et al, 20081 ), field 
lines could easily slip over larger distances than in our simulation, implying that 
the particle beam could sweep over a larger angular range. 

According to |Droge et a/.|(|2010p and lChoUet and Giacalone](|201ip . energetic 
particles should not be scattered significantly during their interplanetary propa- 
gation, implying that energetic particles propagate in the inner heliosphere along 
the same field line. Thus, the association of slipping and interchange reconnection 
provides a new way to understand SEP connection to the Earth and the large 
angular range over which particles can be detected in the interplanetary medium, 
suggesting that the cross-field transport of accelerated particles in the coronal 
magnetic field is not the only way to inject particles in an extended flux tube. 
This result leads us to believe that the MHD model proposed in this study can 
be fully relevant to explain several observations such as the measurement of 
energetic particles over more than 90° of longitude. 

The model that we propose for sweeping particle beams in the heliosphere 
is based on the assumption that the macroscopic spatial evolution of propaga- 
tion channels follows the dynamics of interchanging and slip-running field lines. 
Future work will therefore have to calculate explicitly the particle trajectories 
in the diffusion region and along the continuously reconnecting field lines in a 
time-dependent magnetic field. 
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Figure 1. Cartoon displaying the connectivity problem between the open coronal field lines 
from the acceleration site (green Une) and the Earth-connected parkcr spiral field line (red 
line) . 
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Figure 2. Top left panel: Initial magnetic topology displaying the separatrix field lines in grey. 
At the photospheric boundary, the greyscale shading represents the vertical magnetic field Bz 
at the photosphere, z = 0. The blue field lines forming the dome materialize the fan surface 
and the singular field line intersecting perpendicularly the fan at the null point corresponds to 
the spine. Top right panel : Photospheric distribution of the mesh. The vertical magnetic field 
and the velocity flow in the {x, y) plan. The grey grid displays the non- uniform photospheric 
mesh, the white dashed contours are the negative isocontours of the vertical magnetic fleld 
for Bz = —0.75, —0.5, —0.25, —0.1, at z = 0, and the white solid contours represents the 
positive isocontours of the vertical magnetic field for Bz = 0.5, 0.75, 1 at 2: = 0. The black 
contours displays the photospheric velocity field, orientated toward the negative-y. Bottom 
panel : Vertical two dimensional {x,z) cut of the density distribution \j\, respectively on the 
left and on the right panel, at t = 50 and at the y position of the null point (y = 0). The white 
grid correspond to the mesh in the (x, z) plan at y = 0. (A full-resolution version of this figure 
is available in the online version of the paper) 
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Figure 3. Vertical profiles of the gravity (top left), the temperature (top right), the magnetic 
field (bottom left) and the density (bottom right) at the null point position. 
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Figure 4. (top row) Vertical two dimensionnal (x,z) cut of tiie current density p and of 
the density distribution \j\, respectively on the left and on the right panel, at f = 35 and at 
the y position of the null point (y = 0). (bottom row) The left and the right panel display 
respectively the null point configuration at the initial time and at t = 35. The dark blue 
lines represent the fan and the spine separatrices passing through the null point. The colored 
field lines initially located in the outer connectivity domain are plotted from fixed footpoints 
and the group of colored field lines initially located below the fan surface are plotted from 
footpoints which are fixed in the advected fiow. The distribution of the photospheric vertical 
magnetic field (z = 0) is coded by a greyscale. 
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Figure 5. Evolution of selected field lines. The dark blue field lines are the spine and fan 
field lines passing through the null point. The green field lines represent the open diverging 
magnetic field. The other colored field lines arc plotted with fixed footpoints in the prescribed 
photospheric fiow.The distribution of the vertical photosphcric magnetic field Bz{z = 0) is 
coded in grey scale. Panel (a) : initial [t = 0) magnetic configuration : all colored field 
lines are located below the fan surface in the inner connectivity domain. Panels (b), (c) and 
(d) : evolution of magnetic field lines in the outer connectivity domain. On the first and 
fourth panels, the purple isocontours plotted at the top of the box corresponds to the spatial 
distribution of the logarithm of squashing factor for log Q = 1.12, 4.01 and 6.90. (An animation 
is available in the electronic version of the journal) 
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Figure 6. Evolution of the slipping (solid lines) and plasma (dashed lines) velocity of the 
three selected field lines. The color of the curves corresponds to the color of the selected field 

line. The horizontal line corresponds to the fast-mode speed, cp = + Cg, at the top of 

the box 
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